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! Abstract 

^ ' The linearity of the Hubble relationship(i.e. between m and z) has been tested for 

5t ! galaxies and supernovea for low redshifts. We have studied this relationship for quasars 
for data taken from Veron Cetti Catalogue(2003).The data from Veron Cetti Catalogue 
for quasars appear to be truncated. The data have been analyzed using various statistical 
methods which are suitable for analysing the truncated data. This analysis shows lineraity 
(in logz) of Hubble law for very small z but non-linearity for high redshift.This will shed 
new light not only on the quasar astronomy but also in the cosmological debate. 
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1 Introduction 



The general relationship between a distance of cosmic sources and corresponding redshift 
allows one to establish some important properties for the universe which can be used to 
probe the spatial geometry and especially the underlying cosmological principles. One of 
such relationships is the m-z relationship between the apparent magnitude of the source 
and its redshift. Astronomers normally work with quantities related to apparent magni- 
tude m and absolute magnitude M. The distance modulus is defined as the difference 
between m and M which is related to redshift z and Hq, where Hq is the Hubble con- 
stant. By analyzing the observational data, Hubble formulated a law which states that 
the galaxies appear to be receding with a velocity v proportional to their distance d from 
the observer: 

V — Hod 

. This is known as the Hubble law, where Hq is called the Hubble constant. However, 
this relation can be derived from cosmological theory, if the universe is assumed to be 
homogeneous and isotropic. Various authors^ discussed the limits of validity of the Hubble 
relation. In the above Hubble relation, the distance d should be very large so that the 
recession velocity is larger than the radial component of the peculiar velocities : for 
example this can be up to 1000kms~^ for galaxies inside clusters and this puts a restriction 
, d > 10/i~^Mpc. This means the redshift z has to be much greater than 10~^. However, 
the distance should not be so large that the recession velocity exceeds the speed of light. 
Crudely speaking, one can use the above relation for d « SOOh^^Mpc, or z « 10~^. 
The distance d can be shwon as d ~ ~ 3000/i~^ Mpc for 10~^ > z < 10~^. This 
equation may be considered as first order approximation to the formula for luminosity 
distance as a function of redshift z in the Priedmann model. There are two luminosity 
functions i.e. absolute luminosity and apparent luminosity. Normally astronomers work 
with absolute magnitude M and apparent magnitude m, where d ~ m — M is known 
as distance modulus and is related to the luminosity distance. Therefore, the relation 
between the distance modulus and logz is considered to test Hubble law. The linearity 
of Hubble relationship has been tested for Galaxies and Supernova for low redshifts^. 
Statistical analysis has been done on various samples for galaxies. Sometimes, the samples 
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are constructed subjectively and often samples are taken from Abell's^ catalogue which 
assumes Bubble's law as the selection criterion. Hoessel et al ^ took samples of 116 galaxies 
from the Abell catalogue which supports Hubble law. Segal and NicoU^ took infrared 
astronomical satellite galaxy samples^ and predicted an alternate redshift-distance law 
i.e. z ^ where p = 1,2,3. Some other attempts^ have been made with the IRAS 
1.2Jy Redshift survey. Efron and Petrosian^ studied the viability of various statitical 
tests for truncated data in connection with redshift survey of galaxies and quasars. From 
the plot of redshifts Zi and log luminosities {=yi) for 210 quasars they found the data 
as doubly truncated data. Here, the trunctaion implies that it is not possible to get the 
information regarding the existence of Zi) if it fell outside of the region Ri where, due 
to experimental constraints the distribution of each yi is truncated to a known interval 
Ri depending on Zi. Truncated data may arise in various experiments. McLaren et al^ 
analysed some experimental results in connection with Red Blood Cell volume distribution 
which lead to truncated data. In section II we will describe nonparametric methods to 
find whether or not the apparent magnitudes m arc independent of redshifts Zi for the 
trucated data taken from Vcron Cetti catalogue for quasars. The cosmological implications 
are discussed in section III. 

2 Statistical Analysis of Quasar Data from Veron 
Cetti Catalogue 

One of the main issues in analyzing astronomical data is to answer the following statis- 
tical question. Is a sample of observed points {zi,mi) in the trncated data set of quasar 
survey consistent with the hyopthesis Hq that Zi and are statistically independent ? 
Efron and Petrosian^ investigated this issue in details using a small sub-sample of quasar 
data. Veron Cetti Catalogue^ like other redshift survey provides a pair of measurements 
(zijUii). Various type of observational biases are ignored. One of the most common bi- 
ases is introduced by limiting of the survey. We can write the data set as {zi,mi) 

for i — 1,2,3 n with rrii < mo. The absolute magnitude Mj can be estimated if we 

assume particular cosmological model. The data set {zi, Mi can be reexpressed satisfying 
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atruncated relationship 

Mi < mo - 5 log d + C 

where C is a constant which can be set to zero. Efron et al anlyzed the data from redshift 
survey of 492 galaxies and the magnitude limit tuq = 21.5 of this survey leads to the 
truncated boundary 

Mi < 21.5 -5 log 

The scatter plot of 

mvs log z 

hint that there is truncation in the data. Here, the idea of truncation is used in the sense 
that the observatiopns 

are observable if some condition or mathematical relation is satisfied say, 

log z > am + b 

, for some a&6 . It appears from scatter plot from Veron Cetti data that there exists at 
least one side turncation. By considering a = 3/7 and b = —64/7 we found that there are 
only 18 data points among the 48683 data points for which log 2; > 3/7m — 64/7. So we 
discard these 18 data points and this number is negligible compared to the size of the data 
set and take log 2; < 3/7m — 64/7 as the truncation relationship. In the next step we will 
use the Test of Independence for truncated data as elaborated by Efron and Petrosian. 
Suppose the data consists of a random sample of n pair from the joint distributions of 

data= {xi,yi),i = 1,2, ....n 

For truncated data we assume that pairs {x, y) are observable only if they satisfy the 
truncation relationship 

y < u{x) 

where u{x) is a monotonic function of x. Following Efron and Petrosian we took x — 
m,y — — log z, u{x) = {—3/7)x + 64/7. The test is to accept independence if 

I (data) I < 1.96 
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and we take the rejection probability of the permutation test to be approximately 0.05. 
Here, if we take Wi — IV the test statistic 

|i^(data)| = 704.162 

with 

Wi = 



l^i ^min 

, which leads locally the most powerful test 

|t^(data)| = 875.594 

. In both the cases the extremely large value of ^^^(data) clearly rejects the hypothesis of 
independence. Here, 

p — value ~ 

where 

p — value 

is the maximum level of significance under which null hypothesis (here, independence of 
two variables) is accepted. So z and m are not independednt. 

In the next step, we will try to get best fitting of these data using regression analysis. 
Then we will investigate the conditons under which we can get Hubble relation. The 
scatter plot of z vs. apparent magnitude m is illustrated in Figure l(a,b). Our regression 
analysis shows that we can use the following relation between m and z. 

log(m - 12) = -4.528 + 16.542;^^/^ - 13.891;^^/^ + 3.884;^^/^ 

for z€{0, 7) i.e. for the whole range of z. 

The following observation motivated us to analyse in a different manner. Here, we 
observed that in the region [0.2950; 2.995] , the truncation distribution of the variable 

(m - 12) 
im - 12) 

can be well approximated by beta distribution with parameters a and which are some 
functions of z. Precisely given the truncation, the conditional distribution of 

(m - 12) 
im - 12) 
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can be well approximated by beta(a(2;), (z)) where 

= (7/3) log(^) + 64/3 

.Actually we did the above analysis for 

z = 0.2950, 0.3050, 0.3150, , 2.9950 

. This information has been used to calculate the expected value of m given z for the last 
said values of z. Then wc went on doing usual regression analysis to find out E{m\z) = 
19.484 + 0.886ln{z) - 0.783(log(2))^ for the region [0.2950; 2.9950]. We found the 95% 
tolerance interval with coverage probability 0.95 in the similar fashion. A% tolerance 
interval with coverage probability a means that A% of the future observation will fall in 
the said interval with probability a. In the specified region they are actually {ml{z),mu{z)) 
for given z where 



and 



16.8 + 7.6263^ - 4.162^^ + 0.80^^^ 



mi = 12.51 + 5.576Z - 1.686^^ 



. They are shown in Fig.2. For the region (0; 0.2950) we use our general regression 
techniques to find out our prediction equation as m = 20.060 + 2.139 log 2; and prediction 
interval as 

20.060 + 2.139 log ^1.9631 V0.4573(l. 0122 - 0.1132 logz + 0.2937(log zf 

Here, the prediction interval means that given z the value of m will fall in that interval 
with probability 0.95. 

3 Possible Implications : 

It is found from our analysis for the data from Veron Cetti Catalogue that the Hubble 
relation between m and log;^ is valid for small z i.e. for the range 

z = [0; 0.295] 
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. For higher values oi z, we get different relation as found from regression analysis. 
Conventionally, the Hubble relation is explained as due to the Doopler mechanism for 
shift of the spectral hues. Now the deviation from Hubble relation may be due to some 
other mechanism for redshift. It may be pointed out that the environmental effete for 
the quasars may be taken into consideration to explain thism deviation. This kind of 
environmental effect has been modeled in Doopler like mechanism considered in Dynamic 
Multiple Scattering (DMS)theory^°. This DMS is essntially based on the odea of correla- 
tion indeuced mechanism as discovered by Wolf^^. Finally we have plotted another curve 
in Fig. 3 as 

VeSect = V* =m- MySZ 

. This figure clearly indicates the existence of three different clusters of quasars. It is 
possible to identify these classes of quasars. The detail study of these clusters and its 
implications will be studied in a subsequent paper. 
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Figure Caption 

Figure 1(a). Scatter plot of z vs. m 
Figure 1(b). Scatter plot of m vs. log( 
Figure 2. Regressions for m vs z. 
Figure 3. Scatter Plot of Veff vs z. 



This figure "Figurela.jpg" is available in "jpg" format from: 



http : / / arXi V . org/ p s/ as tro-ph/0605 3 5 6 v 1 



This figure "Figurelb.jpg" is available in "jpg" format from: 



http : / / arXi V . org/ p s/ as tro-ph/0605 3 5 6 v 1 



This figure "Figure2.jpg" is available in "jpg" format from: 



http : / / arXi V . org/ p s/ as tro-ph/0605 3 5 6 v 1 



This figure "Figure3.jpg" is available in "jpg" format from: 



http : / / arXi V . org/ p s/ as tro-ph/0605 3 5 6 v 1 



